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A calculation is made of the maximum steady-stale temperature rise due 
to a small amount of dielectric dissipation in a Luneberg lens which is 
continuously illuminated by a powerful microwave transmitter, when the 
surface of the lens is held at a constant temperature. The temperature 
distribution along the axis of the lens is computed both when the focus is at 
Die surface, and also when it is outside tlie surface at a distance equal to 
one-tenth of the lens radius. The numerical results are given in such a form 
that the maximum temperature rise is easily deduced when the loss tangent 
is any linear function of the refractive index of the lens material. In general 
the maximum steady-slate temperature occurs on the axis at some interior 
point between the center of the lens and the focus. The total power dissipated 
in the lens is also computed. Finally, a brief discussion is given of the time 
scale associated with transient healing of the lens. 

I. INTRODUCTION AND SUMMARY 

When a Luneberg lens is to be used as an antenna for a long-range 
radar, 1 it is important to know how much the lens will be heated by 
dielectric dissipation when it is illuminated by a powerful transmitter. 
This paper presents a calculation of the maximum steady-state tempera- 
ture rise in the interior of the lens, when the surface is held at a constant 
temperature and only a small fraction of the incident power is dissipated. 
Since the maximum temperature rise depends critically on the loss 
tangent of the lens material, and since the loss tangent may vary with 
index of refraction, the results are given in such a form that the maxi- 
mum temperature rise is easily deduced when the loss tangent is any 
linear function of the index of refraction. The index of refraction is 
assumed to vary with radius in a manner appropriate for a Luneberg 
lens of the desired focal length. Numerical computations have been 
made for the case in which the focus is at the surface of the lens, and 
also when it is outside the surface at a distance equal to one-tenth of 
the lens radius. 
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In the idealized problem to be studied, the temperature rise per watt 
of incident power will be the same for a distant transmitter, which 
illuminates the lens with an essentially uniform plane wave, as for a 
transmitter at the focus, with a feed pattern corresponding to a uniform 
plane emergent wave. In either case, if the lens is illuminated from a 
particular direction, and if the index of refraction and the loss tangent 
are functions of the radial coordinate only, the heat source distribution 
and the temperature distribution will be axially symmetric, and in all 
practical cases the maximum temperature will occur somewhere on the 
axis of the lens. Hence to find the maximum temperature we need only 
calculate the temperature distribution along the axis. 

Suppose that the loss tangent of the lens material can be adequately 
represented by a linear function of the index of refraction; thus 

tan 5 = An + B, (1) 

where n is the refractive index and A and B are constants. Then we 
shall show that the axial temperature distribution can be written in the 
form 

TU) = Wk\)[AT A (Z) + BT B (Z)], (2) 

where £ is axial distance in units of the lens radius, with £ = — 1 corre- 
sponding to the plane wave side and £ = +1 to the side nearest the 
focus. P is the total power incident on the lens, k is the thermal con- 
ductivity, and X is the free-space wavelength. The dimensionless func- 
tions T A (%) and T B (£) are given in Table I of Section V and are plotted 
in Figs. 2 and 3 for lenses with normalized focal distances of 1.0 and 1.1, 
measured from the center. Note that the maximum temperature given 
by (2) is independent of the lens radius. 

The foregoing remarks apply to the case in which the lens is illumi- 
nated from a single direction, so that the maximum temperature rise 
occurs on the axis. If the total power P striking the lens comes from 
several different directions, we can deduce upper and lower bounds on 
the maximum temperature rise in the "multiaxial" case from a knowl- 
edge of the temperature distribution along the axis in the "uniaxial" 
case. Since the heat conduction equation is linear, the principle of 
superposition guarantees that the temperature at the center of the lens 
is the same in both cases. Also, the maximum temperature in the multi- 
axial case is less than the maximum temperature in the uniaxial case, 
since the maximum temperature point in the uniaxial case is on the 
axis defined by the incident beam, and this point is not on the other 
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axes in the multiaxial case. It follows that the maximum temperature 
rise in the multiaxial case is at least as great as the rise at the center of 
the lens in the uniaxial case, but not as great as the maximum rise in 
the uniaxial case. The numerical example in Section V indicates that 
in practical situations these two bounds may be so close together that 
a more detailed treatment of the multiaxial case would be superfluous. 
The body of the paper is concerned with the determination of the 
functions T A (%) and 7V£) which occur in (2). In Section II we intro- 
duce a convenient approximation to the index of refraction of the 
Luneberg lens, which is exact if the focus is at the surface, and show 
that under this approximation the ray paths are ellipses. In Section III 
we compute the power flow through every element of a lossless lens, 
and the approximate rate of dissipation of heat, assuming small dissipa- 
tion and a loss tangent of the form (1). An integral representation of 
the temperature along the axis is obtained in Section IV, as well as an 
expression for the total dissipated power. Results of numerical integra- 
tions carried out on an IBM 7090 are given in Section V. Appendix A 
contains a proof that in all practical cases the maximum "uniaxial" 
temperature occurs on the axis, while Appendix B is concerned with the 
nature of the mathematical singularity which occurs in the idealized 
model when the focal point is at the surface of the lens. The time scale 
for transient thermal effects is briefly discussed in Appendix C. 

H. RAY PATHS IN A LUNEBERG LENS 

The path of a typical ray in a TAineberg lens of normalized radius 
unity is shown schematically in Fig. 1. In general the path of a light 




Fig. 1 — Typical ray path in a Luneberg lens. 
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ray in a radially symmetric medium with refractive index n(r) is given 
by 2 



9 ~ d °- ± / ro r[nV -«•!»' 



(3) 



where (r,0) are polar coordinates in the plane of the ray, and k is con- 
stant for a particular ray, being determined, together with the ambiguous 
sign, by the direction of the ray at the initial point (r<j,0o). In the 
special case where the rays are all initially parallel, k is equal to the 
initial distance of the given ray from the axis of the system. 

In principle, (3) may be integrated to find the ray path whenever n 
is a known function of r. For a Luneberg lens with focal point offset 
from the surface, however, the relationship between n and r is given 
by two parametric equations" involving a function defined by a definite 
integral, and it does not seem possible to obtain the equation of the ray 
path explicitly in terms of known functions. An approximation to the 
refractive index which does permit analytic integration of (3) is 

n - [n 2 - (n 2 - l)r 2 ] J , (4) 

where n Q is the index at the center of the lens according to the accurate 
theory; n is a decreasing function of the focal length r . The relation- 
ship (4) is exact if the focus is at the surface (n = V2), and is a good 
approximation if the distance from the focus to the surface is small. 

To find the equation of a typical ray in a lens whose refractive index 
is given by (4), it is convenient first to locate the "turning point" 
(r*,0*) at which the distance of the ray from the center of the lens is a 
minimum. The turning point is defined by 

r n(r ) = k, (5) 

which yields, using (4), 

*/ \ \ n * ~ ["o 4 ~~ 4k 2 ("o 2 - 1)]'"]' 

r w - L - -2<s^Tj J ■ (b) 

The corresponding angle is derived from ( 3 ) , setting 6 = ir when r = «> 
and noting that and r decrease together. We obtain 



w-'+Cra+f 



dr 



r(r* - k 2 )* Jx r[r 2 n 2 (r) - K 2 ]* 

(7) 

Sir . -i , . -l Wo — 2k 

T" sln K ~* tan 2,(1 -««)» ' 
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on substituting (4) into the second integral and carrying out some 
straightforward integrations. 

We can now determine the equation of the ray path inside the lens. 
Proceeding along the ray in either direction from the turning point, we 
have from (3) and (4), 



-»»(«) -*r t ** 

J r - r[nrr 2 — k 2 |* 
= ± \ cos 



2k 2 - rcoV 



(8) 



r 2 lwo 4 - 4k 2 (n 2 - 1)]*' 
or, by rearrangement, 

F(r,d; K ) = r s {[no 4 - 4k 2 (n 2 - l)] 1 cos 2(0 - 0*) + n<?\ - 2k = 0. (9) 

Since 0*(k) is constant for a particular ray, it is easy to write (9) in 
rectangular coordinates and to verify that it is the equation of an ellipse. 

III. RATE OF INTERNAL DISSIPATION OF HEAT 

If the total dissipated power is a small fraction of the incident power, 
as it must be in a practical lens, then the heat losses may be regarded 
as a small perturbation on the power flux in the lossless case, which we 
shall now compute. 

Henceforth we regard each ray as defining a surface of revolution, 
although the ray itself lies in a plane through the axis of the system. 
According to geometrical optics, the total power flow along the tube 
bounded by the ray surfaces corresponding to k and k + <Ik is constant. 
Let dv be the elementary distance normal to the ray in the direction of 
increasing k. If k is regarded as a point function defined by (9), we have 

dK_ BF/B V VF . . 

dv BF/Bk 3F/Bk ' V ' 

where. TF is evaluated by differentiating F with respect to the coordi- 
nate variables while holding k fixed. 

Now let S(r,0) be the power flux along a ray at any point of the lens; 
appropriate units for S with the present normalization of lengths are 
watts/ (radius) 2 . The total power flow along an elementary tube is then 

(IP = 2irr sin 6 S(r,0)dv = constant. (11) 

The constant can be evaluated by considering a ring-shaped element of 
area normal to the incident beam, where k is just the distance of the 
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given ray from the axis. If P is the total power incident on the lens, 
then 

dP = 2P KdK. (12) 

Combining (10), (11), and (12) yields 

flfr ^>- /l«ilK!i - (13) 

The attenuation constant a at any point of the lens is given by the 
well-known approximate relationship 

a = (irn/\) tan 5, (14) 

where n is the refractive index, tan 8 is the loss tangent (assumed 
small), and X is the free-space wavelength, measured here in units of 
the lens radius. Hence the rate of energy dissipation per unit volume, 
in watts/ (radius) , is 

«/ \ « a „\ 2P Kn tan 5 I VF I /lcN 

On the right side of (15), k is defined implicitly by (7) and (9) as a 
function of r and 0. The refractive index n(r) is given by (4), and tan 8 
is supposed to be a known function of n. Differentiation of (7) and (9) 
yields 

— = 2r{W - 4« 2 (« 2 - 1)]* cos 2(6 - 0*) + n 2 }, (16) 

dr 

I B l = -2r{[«o 4 - 4 K 2 (r^o 2 - l)] 1 sin 2(0 - 0*) }, (17) 

r dd 



** = r 2 (-, . ^"^^u cos 2(0 - 0*) 



dK 1 [no 4 - 4ic 2 (no 2 - 1)] J 

+ 2[n 4 - 4, 2 (no 2 - 1)]* sin 2(0 - 0*) S - 4k, 

d0* 1 (W - 2k 2 ) (no 2 - 1) 



(18) 



(19) 



dK (1 - k 2 )* V - 4K 2 (n 2 - 1) 

Hence in principle the rate of heat generation Q(r,6) is a known func- 
tion of position within the lens. 

To determine the power flux from (13) and (16)-(19) when k = 
or k = 1 requires the evaluation of some indeterminate forms. Deriva- 
tions of the following results are straightforward and will be omitted. 
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The power flux along the axis is given by 

S(0 = V L[n*« - (n.» - l)h ~ W - 1) J ' (20) 

where £ is the normalized axial coordinate defined below (2) of Section 
I. The power flux at the equator of the lens is 

lim S(r,ir/2) = P 0/ W. (21) 

r-»l - 

If the focus is at the surface of the lens (n = \/2), then 

lim S(r,0) = 0, < 6 < tt/2. (22) 

r-»l~ 

In this case, particular interest attaches to the power flux in the im- 
mediate neighborhood of the focal point. It is convenient to introduce 
new polar coordinates p,#, with origin at the focus and polar axis in the 
direction of decreasing £ (see Fig. 1). Then for small p the power flux 
is 

S* PoC0 ?*. (23) 

irp 

IV. AXIAL TEMPERATURE DISTRIBUTION AND DISSIPATED POWER 

The steady-state temperature distribution inside the lens satisfies 
Poisson's equation, 

V 2 T = -Q/k, (24) 

where T is the temperature above any convenient reference level, Q 
is the source distribution, and k is the thermal conductivity, expressed 
for the moment in units of watts/ (degree- radius). Since the surface of 
the lens is assumed to be held at a constant temperature, say by air 
conditioning the space between the lens and the radome, the boundary 
condition may be taken as 

T = at r = 1. (25) 

The source distribution Q is a function of the coordinates r,6 only, and 
in all practical cases it decreases with increasing distance from the axis. 
It is proved in Appendix A that the maximum temperature rise then 
occurs on the axis; and as shown in Section I, a knowledge of the axial 
temperature distribution in this case enables us to put upper and lower 
bounds on the maximum temperature rise in a Luneberg lens illumi- 
nated from more than one direction. 
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The temperature distribution on the axis may easily be written down 
in terms of the Green's function for the interior of a sphere; thus, 

T(0 =7 f f Q(r,d)G(r,6^)2Tr 2 sin d dr dd, (26) 

K Jo JO 

where 

G{rM) = (l/4*-)[(r 2 + t - 2r£ cos 0)~ } 

(27) 
- (I + fr 2 - 2r£cos0)~']. 

If desired, one may think of (26) as representing the electrostatic 
potential due to a distributed electric charge of density e Q(r,6)/k 
inside an earthed, conducting sphere of unit radius, but the analogy 
has nothing to do with the mathematics. 

Before starting calculation, we shall assume that tan 5 is a linear 
function of n, say 

tan 5 = An + B. (28) 

Then making use of (15), we may write (26) in the form 

m) = (P /k\)[AT A (Z) + BT B (S)], (29) 

where 



f T f 1 I VF 
T„(S) = 4t / / m -_ 

Jo Jo I dl< I 6k 



G(r,d;Z)r dr dd, 



G(r,0;£)r dr dd. 



(30) 



The dimensionless functions T A (£) and T B (£) are calculated numerically 
in the next section. Note that although the radius of the sphere has 
heretofore been taken as the unit of length, the factor P /k\ has the 
dimensions of temperature, and any consistent set of units (e.g., MKS) 
may be used for P , k, and X in (29). 

It has been tacitly assumed in the foregoing that the thermal con- 
ductivity k is constant throughout the lens. But the conductivity of 
polystyrene foam, out of which Luneberg lenses are usually made, is 
known to increase with increasing temperature, and therefore it may 
be greater in some parts of the lens than in others. However we know 
from general theory that if in a body with a fixed distribution of heat 
sources, the thermal conductivity is increased at any point, the steady- 
state temperature at each point either decreases or remains unchanged. 
Hence the solution of the heat flow problem with a constant value of k 
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less than or equal to the actual value of /." at all points, provides a tem- 
perature distribution which is an upper bound to the actual tempera- 
ture distribution at all points. 

Finally, the fraction of the total incident power which is dissipated in 
the lens is given by integrating Q/Po over the volume. We obtain, from 
(15) and (28), 



AP 
P« 



A AP A + B AP B 



l\ roj 

where a is the radius of the lens in the same units as X, and 






in •'O ^U 



VF 
) | dF/dK 



r dr dd, 

r dr dd. 



(31) 



(32) 



V. NUMERICAL RESULTS 

The functions T A (£) and T B (£) were evaluated on an IBM 7090 by 
a straightforward double application of Simpson's rule to (30). The 
value of k at each point was found by solving (9) by Newton's method; 
then | VF/ (dF/dK) | was calculated from (16)-(19) and G(r,6£) from 
(27). Two values of normalized focal distance, measured from the 
center of the lens, were considered, namely, 

wo = a/2, 



r = 1.0, 
n = l.i, 



n„ = 1.30025. 



(33) 
(34) 



The numerical results are given in Table I, and are plotted in Figs. 2 
and 3. Also the total dissipated power was computed from (32). For 
n = 1.0, 

(AP/Po) = (a/X)[16.9L'l + 13.975]; (35) 

and for 7'n = 1.1, 

(AP/Po) = (c/X)[15.46il + 13.133]. (36) 

For the numerical integration a graded net was used, as follows: 

Region I 



Region II 



r = 0.00 (0.05) 0.40 

6 = 0.0° (7.5°) 180.0° 

/• = 0.40 (0.05) 0.80 

= 0.0° (3.0°) 12.0° 



668 



THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1964 



Region III r = 0.80 (0.02) 1.00 

e = o.o° (i.o°) 12.0° 

Region IV r = 0.40 (0.05) 1.00 
6 = 12.0° (3.0°) 30.0° 

Region V r = 0.40 (0.05) 1.00 
6 = 30.0° (7.5°) 180.0° 

Experimentation with a finer net, obtained by simultaneously halving 
the intervals in r and 6 in Regions III and V, indicates that any errors 
in Table I and (35)-(36) do not exceed a few units of the last figure 
shown. The accuracy is therefore believed to be sufficient for all practical 
purposes. 

A few aspects of the calculation deserve comment. In the first place, 
the Green's function (27) is infinite when the field point coincides with 
the source point; but the integrand of (26) does not become infinite as 
the field point approaches the source point, provided that Q(r,8) re- 
mains finite on the axis. The limiting value of the integrand may be 
either zero or finite, depending upon the direction from which the field 
point approaches the source point; but in any event the contribution of 
the apparent singular point during a naive application of Simpson's 



Table I — The Functions T 


.(£) AND Tb(£) 


« 


ro = 1.0 


ro = 


= l.l 


T A 


T B 


Ta 


Tb 


-1.0 


0.0000 


0.0000 


0.0000 


0.0000 


-0.9 


0.1060 


0.0861 


0.0985 


0.0823 


-0.8 


0.2155 


0.1727 


0.1993 


0.1645 


-0.7 


0.3241 


0.2564 


0.2981 


0.2431 


-0.G 


0.4308 


0.3369 


0.3940 


0.3179 


-0.5 


0.5345 


0.4140 


0.4860 


0.3887 


-0.4 


0.6345 


0.4876 


0.5734 


0.4553 


-0.3 


0.7298 


0.5574 


0.6553 


0.5174 


-0.2 


0.8198 


0.6234 


0.7310 


0.5749 


-0.1 


0.9037 


0.6856 


0.7998 


0.6277 


0.0 


0.9809 


0.7440 


0.8609 


0.6754 


0.1 


1.0493 


0.7974 


0.9123 


0.7169 


0.2 


1.1090 


0.8464 


0.9535 


0.7521 


0.3 


1.1588 


0.8907 


0.9831 


0.7802 


0.4 


1 . 1984 


0.9306 


0.9997 


0.8003 


0.5 


1.2267 


0.9659 


1.0005 


0.8105 


0.6 


1.2386 


0.9937 


0.9785 


0.8051 


0.7 


1.2317 


1.0132 


0.9255 


0.7771 


0.8 


1.2040 


1.0251 


0.8243 


0.7109 


0.9 


1.1521 


1.0324 


0.6206 


0.5552 


1.0 


1.0000 


1.0000 


0.0000 


0.0000 
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Fig. 2 — The functions T ,\ and T B for a lens with focus at the surface (ro = 1). 

rule is at most a finite quantity, which tends to zero as the interval of 
integration is reduced. In the numerical calculation the contribution 
from this point was omitted. 

A discussion of the mathematical singularity at the focus when 
ro = 1 is given in Appendix B. In the mathematical model, the tempera- 
ture near the focus is given by 

T( P ,&) « 7'ocostf, (37) 

where T is a finite constant and p,# are spherical polar coordinates in 
the local coordinate system introduced at the end of Section III. The 
value of r is determined by the source strength in an infinitesimal 
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Fig. 3 — The functions T A and T H for a lens with focus outside the surface 
(r =1.1). 
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region around the focus; in the limit this region makes no contribution 
to the temperature at other points of the sphere. Hence T A {%) and 
T B (£) were calculated, for £ 9^ 1, by setting the integrands of (30) 
equal to zero at r = 1, 8 = 0; and 7^(1) and T B (l) were evaluated as 
in Appendix B. 

It goes without saying that in a physical lens the temperature at the 
focus will be well denned, which the limit of the expression (37) is not, 
and that the temperature distribution will not have an infinite gradient. 
The actual distribution will be determined by the amount of heat that 
the air conditioning can carry away from the immediate neighborhood 
of the focus, as well as by the physical structure of the feed, if the antenna 
is being used for transmission. Since, however, the numerical example 
at the end of this section suggests that the maximum temperature rise 
in a Luneberg lens with surface cooling will be well inside the lens, we 
shall not attempt here a more elaborate analysis of the conditions near 
the focus. 

If ro > 1, so that the focal point is outside the lens surface, then (4) 
is not an exact expression for the index of refraction. However, when 
/•a = 1.1, the maximum difference between the exact index calculated 
according to Ref. 2 and the approximate index is about 0.0055 at about 
r = 0.87, the approximate index being smaller. As a second test, we 
have calculated the distance from the center of the "approximate" 
lens at which various initially parallel rays intersect the axis. The 
distance varies from 1.0881 for paraxial rays (k = 0) to 1.1205 for rays 
with k = 0.95, compared with the design value of 1.1. It tends toward 
infinity for marginal rays, but such rays are insignificant so far as the 
heating problem is concerned anyway, since by hypothesis the surface 
of the lens is in contact with a constant-temperature heat reservoir. We 
therefore feel well justified in using (4) to compute the ray paths for 

r„ = 1.1. 

To give an idea of the size of the numbers involved, Fig. 4 shows plots 
of the axial temperature rise in degrees Fahrenheit per watt of incident 
power, as calculated from (29) and Table I for a lens with the following 
parameters: 

P = 1 watt 

. „ or BTU n nor watts 

A:=0J5 hr-lV-(°F/in) = -° 36 n 1 ^C 

X = 60 cm (500 mc) (38) 

tan 5 = [1 + 25 (n - 1)] X 10 _< 

a = 40 ft 



MICROWAVE HEATING OF A LUNEBERG LENS 



671 



0.07 



H 0.06 
£ 0.05 



Ei 0.04 

I 

z 

111 

cr 0.03 

i 

< 

IL 

ui 0.02 



O 0.01 



-1.0 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1.0 
NORMALIZED AXIAL DISTANCE i 

Fig. 4 — Numerical example of axial temperature distribution in two Luue- 
berg leusea. 
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The assumed values of thermal conductivity and loss tangent are more 
or less representative of polystyrene foam loaded with metal slivers. 
The maximum temperature rise is roughly 0.06°F/watt, and occurs 
about a third of the way from the center of the lens to the surface, in 
the direction of the focal point. It is of interest to note that the maxi- 
mum temperature differs from the temperature at the center of the lens 
by less than 15 per cent, even when the focus is at the surface. 

The power dissipated in the lens is easily calculated from (35) or 
(36). Forr = 1.0, 



(AP/Po) = 0.1777 or 0.85 db; 



(39) 



while for r = 1.1, 



(AP/Po) = 0.1450 or 0.68 db. 



(40) 



Equations (39) and (40) give the dissipated power when a uniform plane 
wave is incident on the lens. Observe, however, that this is not quite 
equal to the power loss when the lens is being used as a transmitter, 
since in that case there is usually a deliberate illumination taper 
across the lens aperture. With a conventional taper, in which the power 
density is higher at the center of the lens than at the edges, the loss 
will be higher than that obtained with uniform illumination; but one 
cannot deduce the total loss from the numbers given in this paper. 
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Note Added in Proof 

A recent paper by Lerner 4 is concerned with essentially the same prob- 
lem as the present paper. On the basis of a considerable number of ap- 
proximations, Lerner calculates the temperature distribution along the 
axis of a surface-focus lens when the thermal conductivity and the loss 
tangent are constant. The "Average" curve of his Fig. 6 is therefore 
comparable to the plot of T B (I-) in our Fig. 2. The two curves do in fact 
lie close together except for points less than a quarter of the lens radius 
distant from the surface focus. Lerner's approximate analysis predicts 
that the maximum temperature rise will occur at the surface focus and 
will be equal to 1.25(P /k\) tan 5. Our calculations give a maximum 
rise of about 1.03(P /kX) tan 5 at a distance of about one-tenth of the 
lens radius from the focus, while the temperature rise at the focus is 
(P /k\) tan 8. 

APPENDIX A 

Position of Temperature Maximum 

We consider the steady-state temperature distribution which satisfies 

V 2 T = -/ for r < a, 

(41) 
T = at r = a. 

The source function / ( = Q/k) is assumed to be axially symmetric 
and nonnegative, and to have continuous first derivatives in the region 
/■ < a. 

We shall use rectangular coordinates (x,y,z), cylindrical coordinates 
(i',<p,z), or spherical coordinates (r,d,<p) as convenient. It is assumed 
that / is independent of <p, and that it is a nonincreasing function of 
distance from the axis, i.e., 

(df/dx) ^ 0. (42) 

Now consider the function 

W = -(dT/dy) = -(dT/dr) sm<p (43) 
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in the hemispherical region S defined by 

r 2 + z ^ a, S V ^ t. (44) 

From (41) and (42) it follows that W is superharmonic in the interior 
of S, since 

V 2 H r = (df/dij) = (df/dv) s\n<p ^ 0. (45) 

Furthermore on the curved surface of the hemisphere we have 

W = -{dT/dy) = -(y/r)(dT/dr) > for ;■ = a, y > 0, (4G) 

since 

(dT/dr) < at r = a (47) 

if / is nonnegative and does not vanish identically. A formal proof of 
this physically obvious statement can be obtained using the Green's 
function representation of the solution of (41). On the base of the 
hemisphere, (43) yields 

W = at // = 0. (48) 

Let U be a harmonic function which takes the same values as W 
on the boundary of S. Since W is superharmonic, 3 we have 

W ^ U (49) 

in the interior of S. But U achieves its minimum value zero only on the 
boundary of S, so neither U nor W can vanish in the interior of S. It 
follows from (43) that dT/dv cannot vanish in the interior of £, and so 
T cannot have a maximum there. Hence the maximum value of T must 
occur on the axis. 

In the present problem the source distribution is given by (15), 
namely 

VF 



f( a\ = Q( r > e) = 2Po tan 5 I Kn 
3K1p) k \k [rsind 



dF/dt 



(50) 



An analytic proof that the right-hand side of (50) is a decreasing func- 
tion of distance from the axis would probably be very laborious. We 
have, however, calculated the expression in braces numerically for the 
two cases treated in this paper, using a square grid of about GOO points 
in r and z, and have verified that on such a grid it is a decreasing func- 
tion of r, except for a very small region near the surface of the lens 
(where r is slightly less than 1 and 6 slightly greater than tt/2 in Fig. 1 ). 
On the other hand, the refractive index n is a decreasing function of r, 
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and in practical cases tan 5 will be a sufficiently rapidly increasing func- 
tion of n so that the whole source distribution will be a decreasing func- 
tion of r. In particular, we have verified numerically that (42) is satisfied 
throughout the lens if tan 8 is given by (38). 

It will be appreciated, of course, that (42) is a sufficient condition but 
not a necessary one to make the maximum temperature rise occur on 
the axis. Whether there could exist a hypothetical loss tangent, having 
a sharp peak at a particular value of ?i, which would lead to a ring- 
shaped temperature maximum instead of to a maximum on the axis is 
still an open question, though not of much practical importance. 

APPENDIX B 

Temperature Distribution near a Surface Focus 

To determine the nature of the temperature singularity in the pres- 
ent mathematical model at a surface focus, we investigate the tempera- 
ture distribution T(p,d) in the half-space # ^ ir/2 due to the source 
function 



Q(p,») = 



C cos t? 



< p < b, 



'- (51) 

0, P ^b. 

Here b and C are constants, and p,d are the polar coordinates introduced 
at the end of Section III. All quantities are independent of the azimuth 
angle (p. We seek a solution which vanishes on the plane d = 7r/2, re- 
mains finite as p — » and as p — * °o , and is continuous, together with its 
normal derivative, at p = b. 
Substitution of a function of the form 

T(p,d) = R(p) cos r? (52) 

into Poisson's equation (24) yields the following equation for R(r): 

*hf)- 2R J- C/k ' °<> <b ' (53) 

dr\ dp/ [ o, p ^ b. 

A solution which satisfies the boundary and continuity conditions is 
easily found to be 



,, -4(l-|), „<„<!,, 



2k 

C_ 
6A-p 2 



R(p) =¥t%, P^h. 



(54) 



MICROWAVE HEATING OF A LUNEBERG LENS 675 

We observe that 

lim R{p) = C/2k, (55) 

and this limit is independent of 6. On the other hand, if p has any fixed 
value greater than zero, 

lim Rip) = 0. (56) 

6-.0 

It follows that the limiting value of the axial temperature is determined 
by the souree distribution in an arbitrarily small hemisphere around 
p = 0, and that in the limit this hemisphere makes no contribution to 
the temperature at other points of the lens. 

If we make the constant C agree with the source distribution given 
by (15) and (23), then using (28) and the fact that n = 1 at the sur- 
face of the lens, we have 

c _ 2aPo _ 2{A + B)P« (5?) 

7T A 

Combining (57) with (55) leads to the results given in Table I, namely 
lim T A (£) = lim T B (Z) = 1. (58) 

APPENDIX C 

Time Scale for Transient Healing 

Suppose that in a Luneberg lens, initially at zero temperature through- 
out, an internal source distribution Q{r,6) is turned on at / = and 
then remains constant in time, while the surface of the lens is held at 
zero temperature. The temperature within the lens must satisfy 

fcV 2 T + Q = pc(dT/dt), (59) 

where k is the thermal conductivity, p the density, and c the specific 
heat of the medium. Conventional units, such as MKS, are used through- 
out this section. 

It is well known that the solution of (59) can be written as the sum 
of a steady-state part and a transient part, i.e., 

T(r,6,t) = T s (r,d) + T t {r,d,t). (60) 

The steady-state term satisfies 

kV 2 T 8 = -Q, (61) 
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while the transient term satisfies 

DV 2 T t = dTt/dt, (62) 

where 

D = k/pc (63) 

is the diffusivity. The transient term vanishes as t — > °o and cancels 
the steady-state solution at t = 0; that is, 

7\(r,0,O) = -T s (r,d). (64) 

At the surface of the sphere both terms vanish; namely, 

T.(a,6) - T t (a,d,t) = 0. (65) 

A transient solution sufficiently general for our needs may be written 
in the form 

00 00 

Tt(r,d,t) = ZZ A nm j n ( Xnm r/a)P n (cos d)e-" nmt , (66) 

„=0 m=l 

where j n is a spherical Bessel function defined in terms of the ordinary 
Bessel function by 

j n (x) = (7t/2.t)'J„ + .(.i-), (67) 

P n is a Legend re polynomial, % nm is the rath root of j„(x), and 

a„m = DXnm /& • (68) 

The Bessel and Legendre functions form complete, orthogonal sets, 
so that at t = any reasonable function of r and 8 may be expanded in 
a double series of the form (66), where the coefficients A nm are given 
by integrals similar to those which define the coefficients in a double 
Fourier series. In particular, if we had calculated the steady-state solu- 
tion T„(r,d) at all points of the sphere, we could expand it in such a 
series and thus satisfy the initial condition (64). We shall not compute 
the coefficients A nm ; we merely observe that in a transient solution of 
the form (66), the individual terms decay exponentially with time, the 
faster the larger a nm . The longest-lived term is the one with smallest 
a, namely 

A Q1 j (xoir/a)e 01 = A 01 / — e K . (69) 

TT/a 
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The rate at which this term decays permits us to define an approximate 
''thermal time constant" for the transient solution, unless of course the 
steady-state solution is such that Aoi is zero or very small compared to 
the other coefficients. Comparing the form of (69) with the expected 
form of the steady-state solution makes it appear obvious that A®, will 
not he unusually small; and therefore an estimate of the time required 
to establish the steady-state solution is furnished by the "half-life" of 
the lowest mode, 

/.,! = 1/aoi = a/rD. (70) 

Assuming for polystyrene foam the numerical values 

p - 1.5 lb/ft 3 , 

c = 0.32cal/gm-°C, (71) 

BTU 

U -°hr.ft 2 -(°F/in)' 

we find after some conversions of units, 

D = 1.12 X 10" 6 m 2 /sec. (72) 

For a sphere of diameter 1 ft, 

/„. = 35 min, (73) 

and for a sphere of diameter 80 ft, 

to, = 156 days. (74) 

The heating time for a large Luneberg lens may thus be several months 
after the transmitter is turned on, with a similar cooling time after it 
is turned off. 
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